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INTRODUCTION 


The  load  interaction  effects  on  fatigue  crack  propagation  due  to  variable 
amplitude  loading  have  been  a subject  of  investigation  for  many  years  and  will 
undoubtedly  continue  to  be  so  for  some  time.  The  inability  to  adequately  pre- 
dict these  interaction  effects  has  prompted  recent  interest  in  studying  the 
interaction  due  to  simple  load  patterns  and  its  influence  on  fatigue  crack 
propagation.  Of  particular  interest  is  the  decrease  in  growth  rate  (delay 
effect  or  crack  retardation)  which  normally  follows  a high  overload  or  a reduc- 
tion in  the  load  level.  This  delay  can  have  a significant  influence  on  the 
fatigue  life  of  a structure  [1,2]^. 

An  understanding  of  the  effect  of  load  interactions,  resulting  from 
changes  in  the  cyclic  load  level,  on  the  fatigue  crack  propagation  is  desirable 
both  in  design  and  in  the  prevention  of  failure  of  components  subjected  to  time 
varying  load  histories.  Several  studies  have  been  undertaken  to  investigate 
the  delay  effect  due  to  simple  load  patterns  such  as  single  and  multiple  over- 
load/underload sequences  [2-18],  Various  explanations  and  models  have  been 
suggested  for  describing  this  delay  behavior  [2-21],  however,  it  is  still  not 
entirely  clear  which  parameters  are  significant. 

Although  reduction  of  the  overload  effect  has  been  observed  for  subsequently 
applied  underloads,  a systematic  study  is  needed  for  the  definition  and  isola- 
tion of  the  associated  variables.  Such  a study  is  essential  before  predictive 
models  for  fatigue  crack  growth  can  be  used  to  adequately  describe  the  influ- 


^ Numbers  in  brackets  refer  to  references  at  end  of  report. 
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ence  of  such  loadings.  This  study  was  undertaken  to  quantitatively  investi- 
gate the  fatigue  crack  delay  or  retardation  effects  of  overstressing  and 
understressing  on  an  otherwise  constant  amplitude  load  cycling  situation.  By 
considering  a three  dimensional  matrix  representation  of  the  selected  experi- 
mental variables,  each  variable  can  be  studied  and  its  influence  on  the  crack 
retardation  phenomenon  assessed.  In  addition  to  isolating  the  effects  of 
the  various  parameters  on  the  delay  behavior,  this  study  was  also  undertaken 
to  determine  if  the  effective  stress  parameters  at  the  crack  tip  could  be 
defined  using  an  inverse  technique.  With  this  technique,  the  growth  rate 
following  the  overload/underload  sequence  is  measured.  This  measured  growth 
rate  is  then  used  in  the  constant  amplitude  growth  rate  equation  to  determine 
an  effective  applied  stress  intensity. 

Using  the  results  of  this  investigation,  an  extended  crack  closure  model 
is  developed  which  accounts  for  all  of  the  quantitative  trends  in  the  data, 
and  also  can  be  used  to  predict  the  crack  opening  stress  intensity  level  and 
the  number  of  delay  cycles  for  a given  set  of  test  conditions.  The  crack 
closure  phenomenon  is  interpreted  as  an  effect,  capable  of  quantifying 
retardation  in  crack  growth,  but  caused  by  the  presence  of  residual  stresses 
in  the  vicinity  of  the  crack  tip. 

Additional  details  of  this  study  are  contained  in  References  22-24. 
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EXPERIMENTAL  INVESTIGATION 


Selection  of  Variables 
and  Development  of  Test  Matrix 


A test  program  was  designed  which  isolated  what  were  believed  to  be  the 
relevant  parameters  affecting  fatigue  crack  growth  in  the  presence  of  single 
overload/underload  cycles.  An  earlier  observation  [11]  indicated  an  apparent 
difference  in  crack  growth  delay  due  to  overloads  which  create  plane  stress 
conditions  and  overloads  which  create  plane  strain  conditions.  Therefore, 
was  chosen  such  that  the  crack  tip  would  be  in  a state  of  plane  stress.  A 
state  of  plane  stress  was  assumed  to  exist  if  the  plastic  zone  diameter,  2r^, 
due  to  the  overload  was  greater  than  the  thickness  of  the  specimen,  B,  where 


-H 

jr‘0  I 
'•  ys-' 
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<1) 


For  this  study,  the  overload  stress  intensity  level  was  chosen  to  be  the  same 


for  all  of  the  tests  ~ 33.33  ksi*^m)  which  gave 

Oiu 


2r  = 0.109  > B = 0.100  in. 

y 


<2) 


Thus,  all  fatigue  cracks  for  the  tests  performed  were  propagated  to  establish 
steady  state  conditions,  overloaded  into  plane  stress  and  returned  to  steady 
state  loading. 

The  parameters  chosen  to  be  investigated  were  tlhe  stress  intensity  ratios 


overload  level  to  maximim  fatigue  level,  underload  level  to 

overload  level,  and  minimum  fatigue  level  to  overload  level, 

R^  = overload/underload  load  sequence  and  the  above  parameters 

are  illustrated  in  Figure  1.  These  parameters  were  believed  to  be  the  signifi 
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cant  loading  parameters  necessary  to  correlate  the  delay  effects.  By  specifying 
the  value  of  each  of  these  ratios  and  one  stress  intensity  level  (in  this  study, 
Kql  the  same  for  all  tests)  the  other  test  conditions, 
are  established. 

These  values  were  selected  because  1)  when  isolated,  would  conveni- 

ently represent  the  effect  accompanying  underloading  from  the  overload  to  the 
underload  value,  a key  parameter  for  both  residual  stress  and  opening  stress 
concepts;  2)  Q had  already  been  shown  to  be  an  important  parameter  in  the 

UJLi 

determination  of  the  arrest/delay  boundary  by  Probst  [8]  and  Himmelein  [11] 
and  3)  had  been  shown  to  be  of  importance  in  tests  with  overloads  only  [11] . 
In  order  that  each  parameter’s  effect  on  fatigue  life  be  properly  assessed,  a 
three-dimensional  test  matrix  was  developed  with  Qql>  forming  the 

coordinate  axes.  Using  such  a scheme,  lines  parallel  to  each  axis  isolated  a 
single  variable  and  an  analysis  of  the  test  results  would  yield  information 
on  how  fatigue  crack  propagation  was  affected  by  that  one  variable.  Two  planes 
parallel  to  the  coordinate  planes  were  chosen  for  principal  investigation: 
vs.  and  vs.  as  shown  in  Figures  2 and  3.  In  addition,  tests 

in  a third  plane  which  spanned  the  values  of  the  and  R^  ratios  used  in 
the  other  two  planes  were  run  to  assess  the  accuracy  of  the  correlations. 

See  Figure  4.  It  should  be  noted  in  Figures  2 and  3 that  these  two  planes 
intersect  at  = 1.8  and  = 0.22  and,  therefore,  the  middle  row  of  tests 
in  each  plane  is  the  same  set  of  tests. 

Since  correlation  of  the  test  results  with  opening  stress  and  residual 
stress  concepts  determined  from  an  inversion  technique  was  desired,  a primary 
objective  in  performing  the  tests  was  to  obtain  measured  crack  growth  rates 
within  the  load  interaction  zone  due  to  the  overload.  Therefore,  precise 
data  of  crack  length  and  number  of  cycles  were  required  following  the  overload. 
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By  choosing  all  of  the  tests  to  have  the  same  overload  level,  the  effect 
of  varying  the  overload  level  was  not  investigated.  The  tests  included  in 
this  study  were  chosen  to  study  the  principle  effects  of  the  selected  para- 
meters. In  all  tests  the  underload  followed  the  overload  with  the  exception 
of  two  tests,  R-12  and  R-13,  (not  shown  in  the  matrix)  which  were  run 
specifically  to  determine  the  effect  of  reversing  the  order  of  the  overload/ 
underload  sequence . 


Test  Specimen  and  Material 

The  test  specimens  used  in  this  study  were  cut  from  2024-T3  aluminum 
alloy,  0.100  inch  thick  from  the  same  stock  as  the  specimens  used  in  Ref.  [11]. 
Using  ASTM  standard  specimens,  the  material  was  determined  to  have  the 
following  tensile  properties  [11]: 

Yield  strength  = 56.9  ksi 
Ultimate  Strength  = 69.0  ksi 
% Elongation  = 14.2% 

The  geometry  for  the  fatigue  test  specimen  is  shown  in  Figure  5.  The  stress 
raiser  shown  in  detail  in  Figure  5 was  machined  with  an  electron-discharge 
technique.  The  specimens  were  obtained  with  a mill  finish  and  were  polished 
to  a mirror  finish  in  the  vicinity  of  the  crack  path  to  facilitate  optical 
observation  of  the  crack  tip  during  crack  growth  measurement.  The  loading 
was  applied  parallel  to  the  direction  of  rolling  of  the  material. 

Test  Equipment 

The  specimens  were  loaded  in  a 20  kip,  closed-loop,  electro-hydraulic 
test  machine  operated  in  load  control.  The  feedback  signal  was  continuously 
monitored  with  a specially  designed  readout  which  provided  a digital  display 
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of  peak  to  peak,  maximum,  minimum  or  mean  level  of  the  cyclic  signal.  A 
counter,  also  connected  to  the  feedback  signal,  was  used  to  count  the  number  of 
cycles  following  the  overload  during  each  test.  This  counter  was  attached 
directly  to  a printer  which  was  used  to  record  the  number  of  cycles.  A strobe 
synchronized  with  the  input  signal  was  used  to  illuminate  the  viewing  surface 
of  the  test  specimen. 

The  crack  length  was  measured  optically  with  a lOOx  microscope  mounted 
on  a two  directional  traversing  system.  A digital  resolver  system  on  the 
horizontal  traverse  produced  a digital  readout  in  the  horizontal  direction 
with  a resolution  of  0.001  mm  (0.00004  in.).  The  direction  of  travel  of  the 
traverse  was  never  changed  during  a test  to  eliminate  any  hysteresis  effects 
in  the  system.  The  data  were  collected  by  advancing  the  microscope  an  incre- 
ment distance  (0.01  or  0.02  mm)  and  recording  the  number  of  cycles  when  the 
crack  had  advanced  that  increment.  For  the  slower  growth  rates,  ^ < 5 x 10  ^ 
in/cycle,  the  ability  to  observe  the  crack  tip  suggests  an  experimental  error 

in  the  measurement  of  Aa  to  be  in  the  order  of  + 0.002  mm  (0.00008  in.).  At 

d 3.  “5 

the  faster  growth  rates,  ~ > 5 x 10  , this  may  be  in  the  order  of  + 0,01  mm 

(0.0004  in.). 

da  “7 

For  very  slow  crack  propagation  rates,  — £ 5 x 10  in/cycle,  the  number 
of  cycles,  N,  was  recorded  every  0.01  mm  (0.0004  in.)  with  all  N values  within 
an  estimated  interval  of  confidence  of  + 50  cycles  of  the  true  value. 

Lightweight  aluminum  compression  guides  lined  with  1/8  in.  felt  were  used 
to  support  the  specimen  when  compressive  underloads  were  applied.  The  guides 
were  removed  for  subsequent  cycling.  The  optical  system  is  shown  in  Figure  6 
and  the  guides,  together  with  a specimen,  are  shown  in  Figure  7. 
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Test  Procedure 


Since  the  scope  of  this  study  strictly  involved  the  effects  of  loading  on 
fatigue  crack  propagation,  care  was  taken  to  control  as  many  other  variables 
as  possible.  All  tests  were  subject  to  nearly  identical  environmental  condi- 
tions of  dessicated  air  and  room  temperature.  All  fatigue  cycling  was  done  at 
20  Hz  except  for  tests  1,  2,  and  3,  which  were  cycled  at  rates  of  10  Hz  due 
to  the  rapid  crack  propagation  for  these  tests.  All  overload  and  underload 
cycles  were  applied  at  0.02  Hz. 

In  order  to  have  a basis  for  comparison  between  tests  in  the  test  matrix, 
the  stress  intensity  factors  were  held  at  quasi-constant  values  throughout 
each  test  by  shedding  the  applied  load.  The  applied  load  was  reduced  at 
every  5%  increase  in  crack  length.  This  insured  that  the  actual  K level  was 
maintained  between  0%  and  +3%  of  the  desired  values.  The  load  levels  were 
calculated  according  to  Tada^s  modification  of  Fedderson's  formula  [25]  for 
the  stress  intensity  factor  of  a center  cracked  specimen.  This  is 


P ^ .K(2b.t)  . _1 

f(a/b) 

K = desired  stress  intensity  factor  (ksii^in) 
t = thickness  of  specimen  (in.) 
a = one-half  crack  length  (in.) 

2b  = width  of  specimen  (in.) 
f(a/b)  *=  correction  factor  for  specimen  geometry 


= [1  - ,025(a/b)^ 


+ .06(a/b)^] 


(3) 


To  make  certain  that  only  single  overload/underload  effects  were  observed, 
care  was  taken  in  each  test  to  establish  equilibrium  before  the  application  of 
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the  overload /underload  sequence.  The  steady  state  condition  prior  to  the 
overload  was  achieved  by  propagating  the  fatigue  crack  at  least  6 iran  (0.24  in.) 
(over  twice  the  theoretical  overload  plastic  zone  diameter)  at  the  quasi- 
constant stress  intensity  values  specified  for  the  test.  The  overload/ 
underload  sequence  determined  from  Equation  3 using  the  measured  crack  length 
was  then  applied.  Following  this,  cyclic  loading  was  resumed  at  the  pre- 
overload level  with  the  load  being  shed  in  the  same  manner  as  described  above 
to  maintain  the  quasi-constant  loading.  The  crack  length,  a,  at  the  point  of 
the  overload  application  varied  between  0.27  mm  and  2.03  mm  (0.41  to  2.03  in.). 
The  corresponding  overload  net  section  stress  varied  between  0.27  and  0.55  of 
the  yield  strength  of  the  material. 

Cycling  was  continuous  for  both  pre-overload  and  post-overload  fatigue 
cycling  to  negate  the  possibility  of  time  or  underloading  to  zero  load  affect- 
ing subsequent  crack  growth. 

The  crack  length,  a,  and  number  of  cycles,  N,  were  monitored  continuously 
for  each  test  and  discrete  data  points  collected  every  0.01  mm,  0.02  mm,  or 
0.05  mm  (0.0004,  0.0008  or  0.002  in.),  the  larger  intervals  used  for  more 
rapid  crack  propagation  rates.  Data  were  recorded  by  advancing  the  optical 
system  by  the  specified  increment  and  pressing  a trigger  button  to  the  printer 
(monitoring  number  of  cycles)  when  the  crack  had  grown  to  the  incremented 
position.  Typical  fracture  surfaces  are  shown  for  three  specimens  in  Figure  8. 
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DATA  REDUCTION  - NUMERICAL  DIFFERENTIATION 


One  of  the  objectives  of  this  study  was  to  characterize  the  growth  rate 
behavior  through  the  overload  affected  zone  following  the  overload /underload 
sequence.  To  obtain  the  growth  rate,  da/dN,  it  was  necessary  to  record  suf- 
ficient a versus  N data  and  then  differentiate  these  data.  Several  different 
techniques  have  been  used  for  the  numerical  differentiation  of  experimental 
data.  Incremental  differentiation  methods  tend  to  amplify  any  small  scale 
variation  in  the  actual  data.  As  an  alternative,  a spline  function  method 
was  selected  in  which  a series  of  cubic  polynomials  were  fit  to  the  experi- 
mental data  and  then  analytically  differentiated  to  obtain  the  crack  growth 
rate.  This  method,  frequently  used  to  differentiate  experimental  data  [26,27], 
provides  a smoothed  crack  growth  rate  curve. 

To  apply  the  spline  method,  third  order  pol3momials  were  fit  to  three 
intervals  of  the  entire  post-overload  a versus  N data  of  a given  test.  The 
values  of  the  polynomials  and  their  first  two  derivatives  were  matched  at  the 
knots  (points  where  the  polynomials  join).  The  computer  technique  of  de  Boor 
and  Rice  [28,29]  was  modified  to  allow  the  first  interior  knot  to  be  fixed  at 
a specified  position.  This  first  knot  was  located  after  the  intlal  decelera- 
tion where  the  a versus  N curve  started  to  level  out.  It  was  found  for  all 
tests  that  this  knot  could  be  located  at  a point  equal  to  5 per  cent  of  the 
calculated  plastic  zone  diameter.  The  computer  routine  then  optimally  located 
the  second  knot  and  determined  the  coefficients  of  the  three  polynomials  which 
minimized  the  least  square  error  over  the  entire  range  of  the  data. 

Figure  9 shows  a typical  fit  of  the  spline  function  to  the  a versus  N 
data  and  Figure  10  shows  the  corresponding  analytically  differentiated  curve, 
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i.e.,  — versus  a.  The  spline  method  of  differentiation  was  compared  to  a 
linear,  seven  point  least  squares  movable  strip  method  [2,30],  The  resulting 
growth  rate  from  this  method  is  also  shown  in  Figure  10  for  the  same  data. 

To  compare  the  two  methods,  the  da/dN  from  the  movable  strip  method  was  numeri- 
cally integrated  and  compared  with  the  original  data.  This  is  shown  in  Figure 
9.  With  the  spline  method,  integration  of  da/dN  will  simply  reproduce  the 
original  spline  function  since  the  derivative  is  obtained  analytically.  As  can 
be  seen  from  Figure  9 and  10,  the  spline  method  provides  a good  representation 
of  the  a versus  N data  and  a smoothed  da/dN  curve.  An  analytical  expression 
is  also  obtained  for  an  entire  data  set.  It  should  be  pointed  out  that  with 
the  movable  strip  technique  a considerably  smoother  differentiation  than  that 
shown  can  be  obtained  by  using  a wider  strip  (either  by  increasing  the  number 
of  data  points  or  by  using  only  every  third  or  fourth  data  point) . 

One  of  the  important  characteristics  of  the  delay  behavior  is  the  minimum 
growth  rate  following  the  overload.  Therefore,  the  minimum  growth  rate  for 
each  test  as  determined  from  the  spline  method  was  compared  with  the  data 
graphically.  In  some  of  the  tests  the  spline  method  produced  a slightly  lower 
value  for  the  minimum  growth  rate.  For  these  tests  (numbers  3,  10,  17,  23, 

32,  33,  34)  the  graphical  values  were  used. 
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TEST  RESULTS 


Constant  Amplitude  Tests 


For  the  purpose  of  comparing  the  fatigue  crack  growth  rates  obtained  from 


the  actual  test  data,  two  constant  stress  amplitude  tests  were  performed  at 


each  of  five  different  fatigue  stress  ratios:  R^  = 0.01,  0.2,  0.33,  0.44, 

0.55.  By  performing  constant  stress  amplitude  tests,  the  value  of  the  stress 
intensity  factors  increased  with  increasing  crack  length,  while 

Rp  was  constant.  By  measuring  crack  length  and  the  niimber  of  cycles,  the 
crack  growth  rate  could  be  determined  for  a given  value  of  AK  and  corresponding 


Rp.  The  results  are  shown  in  Figure  11. 

A series  of  three  cubic  splines  [28,29]  were  fit  to  each  set  of  a vs.  N 
(approximately  85  points  per  set)  with  the  two  interior  knot  positions 
optimized  so  as  to  minimize  the  least  square  error.  The  values  of  da/dtl  were 
obtained  by  differentiating  the  resulting  spline  functions.  Mth  a knowledge 
of  the  loads,  ^MIN*  well  as  the  crack  length  and  number  of  cycles 

as  represented  by  the  spline  function,  the  values  of  AK^pp  and  the  corresponding 
growth  rate  were  determined  at  170  points  for  each  test,  yielding:  a total  of 
1700  points  for  all  ten  tests.  From  these  values  of  da/dN  and  AK  , the 
values  of  the  two  constants  for  the  steady  state  growth  rate  equations  were 
determined  using  a least  square  fit.  This  gave  for  Elher's  equation  |19,20| 


i = 

AKppp  = (0.5  f 0.4Rp)AK 

C = 6.92  X 10“^ 
n = 3.86 
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The  data  were  also  shown  to  agree  with  Gallagher’s  equation  [31] 


||  = A[AK(1  + 3Rp)]^  (5) 

where  3 = 0.698  for  R > 0 

r 

3 = 0.418  for  R < 0 

r 

A = 6.78  X 
b = 3.58 

The  constants  in  this  equation  were  determined  by  Gallagher  [31]  using  Hudson's 
data  [32]. 1 

Equation  4 was  used  in  determining  the  experimental  parameters  associated 
with  extended  crack  closure  and  Equation  5 was  used  in  calculating  the  residual 
stress  intensity  parameters.  These  concepts  are  discussed  in  detail  later. 

Overload /Underload  Tests 

A summary  of  the  test  conditions  and  results  for  all  of  the  overload/ 
underload  tests  is  presented  in  Table  I.  The  actual  loading  conditions  for 
each  of  the  tests  in  the  three  planes  of  the  original  test  matrix  are  shown  in 
Figures  2-4. 

Following  the  completion  of  the  test  program  as  presented,  three  addi- 
tional tests  were  run  in  order  to  better  define  crack  growth  trends  for  tests 
having  negative  R^^  and  high  values , which  were  the  two  parameters  found 
to  be  of  most  significance  in  determining  the  amount  of  delay.  These  tests 
(16-2.4,  17-2.4,  and  16-2.6)  are  in  the  same  vs.  R^^  plane  as  the  original 
matrix  tests  shown  in  Figure  2.  To  better  define  an  arrest /delay  boundary  for 

^ Using  Gallagher’s  form  of  the  growth  equation  a minimum  least  squared  error 
was  obtained  for  the  data  of  this  study  with  A = 6.28  x 10“10,  b = 3.79,  and 
3 = 0.647  (for  R^  > 0 only). 
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negative  values,  three  additional  tests  were  run  with  - 3.0* 

The  data  obtained  from  the  tests  in  the  matrices  illustrate  the  same 
qualitative  trends  as  noted  by  other  investigators.  By  representing  the  a vs. 
N data  from  each  of  the  tests  with  the  spline  functions,  the  results  can  be 
illustrated  by  plotting  the  resulting  spline  functions.  These  results  are 
shown  in  matrix  form  in  Figures  12  and  13  for  the  two  principle  planes. 
Analytical  differentiation  of  the  spline  functions  yields  the  corresponding 
da/dN  vs.  a curves  shown  in  Figures  14  and  15. 

The  effects  of  varying  each  of  the  selected  parameters  (Qqj^>  ^0* 

K^)  with  all  other  parameters  held  constant  can  easily  be  compared  in  these 
figures.  Horizontal  directions  from  left  to  right  on  each  of  the  matrices 
correspond  to  an  increase  in  from  a value  of  -1.00  to  a value  at  which 
there  is  no  applied  underload.  From  the  results,  it  is  seen  that  as  R^^  is 
increased,  that  is,  as  the  underload  value  approaches  the  minimum  stress 
intensity  value,  an  increasing  amount  of  crack  delay  results.  Furthermore, 
it  may  be  seen  from  the  two  matrices  showing  the  growth  rates  of  each  test 
(Figures  14  and  15)  that  as  R^^  increases,  the  minimum  growth  rate  decreases. 

, for  each 


The  number  of  delay  cycles,  and  the  minimum  growth  rate, 


MIN 


test  are  shown  in  matrix  form  in  Figures  16  and  17.  Consideration  of  vertical 
lines  in  Figure  16  yields  information  accompanying  the  variation  of  the  para- 
meter Q for  all  other  parameters  held  constant.  The  five  vertical  lines 
illustrate  that  as  increases,  or  as  the  value  of  decreases  for  other 

stress  intensity  levels  fixed,  the  number  of  delay  cycles  Increases,  with  a 
corresponding  decrease  in  the  minimum  growth  rate.  Examination  of  the 
vertical  lines  of  the  vs.  R^^  matrix  (Figure  17)  does  not  illustrate  as 
significant  a trend  due  to  varying  as  is  observed  for  varying  R^^  or 
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These  observations  suggest  that  is  of  secondary  importance  in  the  predict 
tion  of  delay  as  compared  with  the  significant  effect  due  to  and  Qqj^- 

Efforts  to  correlate  the  measured  overload  affected  zone,  Aa*,  with  the 
selected  set  of  parameters  did  not  yield  any  consistent  trends.  However,  the 
average  value  of  Aa*  was  0.104  in.  (with  a standard  deviation  of  0.029) 
which  agrees  reasonably  well  with  the  overload  plastic  zone  diameter  of  0.109 
in.  determined  from  Equation  1.  This  agreement  between  the  calculated  over- 
load plastic  zone  size  and  the  test  data  is  in  contradiction  to  the  findings 
of  Lankford  and  Davidson  [14]  who  showed  through  the  use  of  electron  channeling 
patterns  that  fatigue  crack  retardation  is  not  related  to  the  maximum  plastic 
zone  dimension.  However,  it  should  be  mentioned  that  since  K was  the  same 

ULi 

for  all  tests,  a range  of  overload  plastic  zone  sizes  was  not  considered  in 
this  study.  The  actual  values  of  the  measured  delay  zones,  Aa*,  were  determined 
from  enlarged  plots  of  the  a vs.  N data.  The  ratio  of  Aa*  to  the  plane  stress 
plastic  zone  size,  2r^,  (2r^  = 0.109  in.)  is  shown  for  each  of  the  original 

matrix  tests  (Figures  2 and  3)  in  Figure  18. 

Tests  R-12  and  R-13  were  run  with  the  identical  test  conditions  as  tests 
12  and  13  except  the  underload  preceded  the  overload.  for  these  two 
reversed  tests  were  160,000  and  169,000  cycles  respectively  as  compared  to 
14,400  and  42,600  cycles  for  tests  12  and  13  respectively.  Test  15  had  the 
same  load  levels  except  with  no  underload  applied  *SiAX  ^IN 

the  same).  for  this  test  was  142,000  cycles  which  compares  favorably  with 

for  the  two  reversed  tests.  This  comparison  indicates  that  an  underload 

preceding  an  overload  has  a different  effect  than  when  the  underload  follows 
the  overload.  Also,  it  suggests  that  if  an  overload  follows  the  underload, 
the  overload  essentially  eliminates  any  effect  due  to  the  underload.  Note 
that  this  is  based  only  on  the  results  of  two  tests  and  should  be  examined  further. 
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INVERSE  METHOD  OF  DETERMINING 
CRACK  TIP  STRESS  INTENSITY  PARAMETERS 

James  and  Anderson  [33]  used  fatigue  crack  growth  rate  data  and  an 
inverse  method  to  determine  the  stress  intensity  factor  for  a complex  geo- 
metry. For  their  method,  the  constant  amplitude  growth  rate  was  determined  as 
a function  of  the  applied  K for  a given  material  using  a geometry  for  which  the 
stress  intensity  factor  was  known.  The  fatigue  crack  growth  rate  in  the  com- 
plex specimen,  made  of  the  same  material,  was  then  measured.  The  applied  K 
which  produced  this  same  growth  rate  in  the  known  specimen  is  then  the  stress 
intensity  factor  (at  the  measured  crack  length)  for  the  complex  geometry. 

A similar  inverse  approach  was  used  in  this  study  to  evaluate  the  change 
in  the  effective  stress  Intensity  parameters  at  the  crack  tip  resulting  from 
the  overload/underload  sequence.  Furthermore,  it  was  possible  to  determine 
how  these  parameters  changed  throughout  the  overload  affected  zone.  The 
inverse  method  Involved  experimentally  determining  the  fatigue  crack  growth 
rate  through  the  overload  affected  zone  and  then  determining  the  effective 
applied  stress  intensity  level  from  the  constant  amplitude  growth  rate  equa- 
tion. This  is  described  further  below. 


Residual  Stress  Intensity  Parameters 


The  minimum  value  of  an  effective  residual  stress  intensity  parameter 
was  detemined  for  each  of  the  tests  as  described  below.  Using  the  following 
form  of  the  constant  amplitude  growth  rate  equation 


^ = A[(l  + SR)AK]^ 


(5) 
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and  solving  for  R gives 


R 


= i /JL  fi  • il 

3 \aK  [a  dNj  " J 


(6) 


The  effect  of  a residual  stress  ahead  of  the  crack  created  by  the 

overload/ underload  sequence  would  be  to  change  the  maximum  and  minimum  values 

of  the  applied  stress  intensity  levels  by  the  amount  of  the  effective  residual 

stress,  [31].  If  the  residual  stress  is  compressive,  then  the  effective 

value  of  the  stress  ratio,  R™„  would  be 

hr  r 


= ^IN  ""R 

^ Vx  - \ 


- K, 


(7) 


Solving  for  gives 
K 


’SlIN  SlAX 


1 - R 


(8) 


EFF 


At  the  minimum  growth  rate  following  the  overload/underload  sequence, 

MAX 

Rgpp  will  be  a minimum.  Therefore,  becomes 

.MIN 


^MAX  _ ^IN 


R 


EFF 


SlAX 


1 - 

^FF 


(9) 


MIN 


At  this  minimum  growth  rate,  be  determined  from  Equation  6 and  the 

experimental  results.  This  gives 


MIN  ^ 1 / j,  fl  . da  lb  _ 1 

^FF  3 \aK  [a  dN  MInJ  ^ 


(10) 
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where 


da 

dN 


MIN 


is  the  minimum  growth  rate  determined  from  the  test  results  and 


MIN 


(is. 

AK  is  the  applied  stress  intensity  range  following  the  overload.  Using  ^ 

MIN 

from  the  experimental  results  as  given  in  Table  I,  was  determined  for 

iir  t 

each  of  the  tests  using  Equation  10.  The  corresponding  values  for  were 

K 

determined  using  Equation  8 where  K is  a compressive  residual  stress.  These 
results  are  presented  in  Table  II. 

It  should  be  noted  that  Table  II  varies  from  -1.51  to  +0.65  with 

most  of  the  values  being  negative.  These  values  were  obtained  using  the  in- 
verse method  and  the  constant  amplitude  growth  rate  equation.  Equation  5, 
however,  the  value  of  3 for  negative  R ratios  in  this  equation  is  based  on  a 
limited  amount  of  data. 

It  has  been  shown  that  the  significant  parameters  associated  with  the 
Wheeler  and  Willenborg  models  are  the  parameters,  and 

respectively  [13,18].  Attempts  to  correlate  the  residual  stress  intensity 
parameter  with  these  parameters  yielded  no  apparent  correlation.  See 

Figures  19  and  20.  This  was  anticipated  since  neither  of  these  models  account 
for  the  effect  of  changing  the  level  of  the  underload.  In  this  study,  a 
significant  effect  resulted  from  changing  the  underload  level. 


Crack  Opening  Stress  Intensity  Parameters 

According  to  Elber’s  crack  closure  theory  [19,20],  a fatigue  crack  remains 

closed  during  a portion  of  the  tensile  part  of  the  load  cycle.  The  effective 

stress  intensity  range,  AK  , that  causes  crack  propagation  is  then 

Er  r 


AK. 


EFF 


^SlAX 


- K 


OP 


(11) 


where  K^p  is  the  stress  intensity  level  at  which  the  crack  opens.  Elber*s 
crack  growth  equation  is  then 
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(4) 


da 

dN 


Substituting  Equation  11  into  Equation  4 and  solving  for  K^p  gives 


^■OP  ^SlAX 


r 

1 

C 


1 

da|n 

dN 


(12) 


Using  the  inverse  method  described  above,  K^p  following  the  overload/ 
underload  sequence  can  be  found,  since  Is  constant  after  the  overload/ 

underload  and  da/dN  was  measured.  da/dN  through  the  overload  affected  zone 
was  determined  from  the  spline  function  representation  of  the  a versus  N data. 
Using  this,  K^p  was  determined  through  the  overload  affected  zone  for  each  of 
the  tests  in  the  two  principle  planes  of  the  test  matrix.  These  results  are 
shown  in  Figures  21  and  22. 

Also,  the  maximuni  value  of  K^p  was  determined  from 


j^MAX  ^ 


SlAX 


fl 

da 

Ic 

dN 

MIN 

J 

(13) 


where  the  minimum  growth  rate  used  in  this  equation  is  the  value  given  in 
Table  I.  These  results  are  given  in  Table  III. 
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EXTENSION  OF  CRACK  CLOSURE 


Elber^s  crack  closure  concept  [19,20]  has  been  used  to  qualitatively 
explain  the  delay  behavior  following  an  overload  condition  [5,21].  In  com- 
paring the  effective  opening  stress  intensity  factor  following  the  overload/ 
underload  sequence  with  the  constant  amplitude  behavior  it  is  assumed  that  the 
overload/underload  sequence  establishes  the  residual  stress  field  in  the 
vicinity  of  the  crack  tip  and  that  subsequent  fatigue  cycling  does  not  change 
this  stress  field.  Also,  it  is  hypothesized  that  the  resulting  crack  closure 
behavior  is  due  to  this  residual  stress  field. 

From  crack  closure,  the  effective  stress  intensity  factor  is 

^^EFF  ^SlAX  " 

and 

§ - '=<‘SlAX  - (W) 


For  constant  amplitude  loading.  Fiber  defined 
U = ~ ^OP 

SlAX  “ SlIN 

and  found  for  2024-T3  aluminum  alloy 


(15) 


U = 0.5  + 0.4  Rp  (16) 

From  the  assumption  that  the  overload/underload  sequence  establishes  the 
residual  stress  field,  U for  the  overload /underload  sequence  can  be  defined 
analogous  to  Equation  15  as 
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(17) 


U 


_ OL  ^QP 
" ^OL  - ^JL 


As  described  in  the  previous  section  was  determined  using  the 

experimentally  determined  minimum  growth  rate  and,  therefore,  U given  by 
this  equation  also  corresponds  to  the  value  at  the  minimum  growth  rate.  If 
crack  closure,  as  presented  for  constant  amplitude  loading,  describes  the 
results,  then  analogous  to  Equation  16 


“oL  ' + “•'"Vo 


(18) 


should  agree  with  the  experimentally  determined  values  for  U . This  is 

UJj 

compared  in  Figure  23  which  shows  that  Equation  18  does  not  adequately 
describe  the  results.  However,  it  is  apparent  that  U ^ is  related  to  the 
loading  conditions,  i.e., 

"oL  ■ %L-  ■ ■ •> 

Examination  of  the  data  suggested  that  this  function  should  be  of  the  form, 

"oL  ' “l  + “z’Vo  + + “4V  <2“) 

Using  the  results  from  the  tests  that  did  not  arrest,  a multiple  linear 
regression  analysis  was  performed  to  determine  the  constants  of  Equation  20. 

= 0.408 
C = 0.367 

(21) 

= 0.117 

C,  = 0.075 
4 
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This  equation  can  be  compared  with  the  experimental  results  by  plotting 

(U  - C,Q  ) versus  This  is  seen  to  agree  very  well  with  the  test  results 

Uti  ^ UJu  uu 

in  Figure  24. 

The  influence  of  was  also  examined  by  writing  Equation  20  as 


“OL  - + <^2^0  + SV  <=4^  ^ 


(22) 


and  performing  the  regression  analysis.  Although  this  provided  a slight 
improvement  in  the  results,  an  analysis  of  variance  showed  that  the 
parameter  when  included  in  this  equation  was  not  statistically  significant. 
A similar  analysis  eliminating  the  term  in  Equation  20  showed  that 
was  statistically  significant. 

By  combining  Equation  17  with  Equation  20  and  solving  for  gives 


+ (23) 

The  values  for  calculated  from  this  equation  are  compared  with  the 

experimentally  determined  values  in  Figure  25.  As  can  be  seen,  the  maximum 

opening  stress  intensity  value  is  predicted  almost  exclusively  within  5% 

MAX 

of  the  measured  value.  Calculating  K^p  from  Equation  23  and  using  the  growth 
rate  equation,  Equation  14,  the  minimum  growth  rates  following  the  overload/ 
underload  sequences  were  determined.  These  results  are  compared  with  the 
experimental  results  in  Figure  26.  The  calculated  and  experimental  results 
for  the  extended  crack  closure  parameters  are  presented  in  Table  III. 

The  effects  of  the  test  parameters  on  the  minimum  growth  rate  can  be  more 
readily  illustrated  by  substituting  from  Equation  23  into  Equation  14  and 


21 


plotting  the  results  from  constant  values  of  Figure  27  shows  several 

d.  a 

level  curves  of  constant  ~ as  functions  of  and  In  this  study 

if  the  growth  was  less  than  0,0004  in.  in  10  cycles,  the  test  was  stopped 
and  the  crack  was  assumed  to  have  arrested.  Therefore,  a minimum  growth  rate 
of  4 X 10  in/cycle,  which  is  the  experimental  resolution  from  the  data, 
can  be  considered  arrest.  This  is  the  curve  labeled  arrest  in  Figure  27. 

This  arrest/delay  boundary  was  based  on  the  finite  growth  test  results,  however, 
it  compares  favorably  with  the  experimentally  determined  arrest /delay  boundary 
found  by  Himmelein  [11]  which  was  for  overloads  with  no  underload  (R^q  ~ Rj^)  • 

The  arrest  tests  from  this  study  are  shown  in  Figure  27.  The  arrest/delay 
boundary  from  Fiber’s  crack  closure  equation  is  also  shown  in  this  figure. 

In  evaluating  the  validity  of  the  extended  crack  closure  results,  it  is 
important  to  keep  in  mind  the  number  of  tests  and  the  distribution  of  the 
loading  conditions  of  these  tests.  Of  the  twenty- four  tests  used  in  deter- 
mining the  coefficients  of  Equation  20,  only  three  tests  had  Q = 2.2. 

UXj 

Furthermore,  eleven  of  the  tests  were  performed  in  the  small  region  of 
0 RyQ  j<  0.3  and  1.6  £ ;<  1.8.  Even  with  this  very  unequal  distribution 

of  the  data,  the  predicted  results  agree  reasonably  well  with  the  actual  data. 
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CORRELATION  OF  THE  NUMBER 


OF  DELAY  CYCLES 

There  are  a large  number  of  parameters  that  Influence  the  delay  behavior 
resulting  from  a change  in  load  level  or  load  sequence.  In  this  study 
only  the  simple  overload/underload  sequence  was  investigated.  For  these 
results  it  was  shown  that  the  maximum  opening  stress  intensity  can  be  corre- 
lated with  the  two  parameters  and  This  extension  of  crack  closure 

can  also  be  used  to  correlate  the  number  of  delay  cycles  with  the  observed 
results. 

From  the  experimental  results  it  was  observed  that  the  number  of  delay 

dal 


cycles,  Nj^,  was  directly  related  to  the  minimum  growth  rate, 
correlation  is  shown  in  Figure  28.  From  Equation  13 


MIN* 


This 


da 

dN 


MIN 


and  the  correlation  of  Figure  28  indicated  that  an  average  growth  rate, 


(24) 

2r 


N. 


would  correlate  with 


*SlAX 


K 


MAX 
OP  ’ 


D 


x.e. , 


2r 


(25) 


where 


2r  = - 
y TT 


a 

V y . 


,,MAX 


MIN 

EFF 


A'  and  b ' were  found  using  a least  squares  fit  on  logarithmic  coordinates , 
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from  Equation  23 


A’  = 8.17  X 10 
b’  = 2.82 

MA 

Solving  Equation  25  for  and  substituting  in  for  K^p 

gives 


N 


D 


[1  - (1  - (C^  + + CjRJ^ 


■s 

» 


(26) 


The  values  of  calculated  from  this  equation  are  compared  to  the  experi- 
mental values  in  Figure  29.  As  can  be  seen  this  shows  good  correlation. 

It  should  be  noted  that  2r^  was  the  same  for  all  tests  in  this  investiga- 
tion, however,  Probst  [8]  also  showed  that  a similar  effective  stress  parameter 
correlated  with  the  average  growth  rate. 

It  is  interesting  to  note  that  the  influence  of  the  load  level  can  be  seen 
from  Equation  26.  By  substituting  in  the  expression  for  r^  and  rewriting  this 
equation  gives 


N. 


.2-b’ 


D 


^OL  ^0’ 


-b’ 


(27) 


Since  f(Q^P^,  is  a function  of  non-dimensional  parameters,  the  load 

level  dependency  is  described  entirely  by  the  term  in  Equation  27*  The 
value  of  b’  is  2.82  and,  therefore,  the  number  of  delay  cycles  is  proportional 
to  raised  to  the  -0.82  power,  i.e.. 


-.82 

OL 


(28) 


This  indicates  a minor  dependency  on  the  magnitude  of  the  load  level  and  that 
the  ratios  of  the  load  levels  are  more  significant.  This  substantiates  the 
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work  of  Himmelein  [11]  who  randomized  load  levels,  but  systematically  varied 
the  load  ratios. 
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CONCLUSIONS 


From  the  results  of  this  study  the  following  conclusions  are  made: 

1)  Utilizing  a three  dimensional  test  matrix,  the  delay  effects  due  to 

each  of  the  non-dimensional  load  parameters,  Qq^>  \[0’  ^ 

could  be  observed. 

2)  By  carefully  measuring  the  crack  growth  rate  through  the  overload 
affected  zone,  the  inverse  method  can  be  used  to  determine  the  crack 
tip  stress  parameters.  This  provides  an  effective  method  for 
studying  the  delay  behavior. 

3)  Utilizing  the  crack  tip  stress  parameters,  it  was  found  that  the 
loading  parameters  and  had  a significant  effect  on  the  delay 
behavior,  the  magnitude  of  the  load  pattern  was  of  secondary 
importance  and  the  influence  of  was  insignificant. 

4)  Comparing  the  results  of  two  reversed  tests  in  which  the  underload 
preceded  the  overload  showed  that  the  subsequent  overload  apparently 
eliminates  any  change  in  growth  rate  behavior  due  to  the  underload. 

5)  The  number  of  delay  cycles  was  found  to  correlate  directly  with  the 
minimum  growth  rate  following  the  overload /underload  sequence. 

6)  An  extension  of  the  crack  closure  concept  was  developed  which 
quantitatively  describes  the  delay  behavior  following  the  overload/ 
underload  sequence.  From  this  the  maximum  value  of  the  opening 
stress  intensity  and  the  minimum  growth  rate  following  the 

over load /under load  can  be  determined.  The  number  of  delay  cycles 
can  also  be  predicted. 
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TABLE  I 


TEST  RESULTS 


Test 

*SlAX 

AK 

''f 

n 

No.  (ksi/TF) 

(ksl  /Tin) 

^OL 

1 

33.33 

11.67 

0.65 

1.00 

2 

25.64 

11.50 

0.55 

1.30 

3 

20.83 

11,50 

0.45 

1.60 

4 

17.54 

11.50 

0.34 

1.90 

51 

15.15 

11.50 

0.24 

2.20 

6 

20.83 

13.50 

0.35 

1.60 

7 

20.83 

13.50 

0.35 

1.60 

8 

20.83 

13.50 

0.35 

1.60 

9 

20.83 

13.50 

0.35 

1.60 

10 
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1.60 

11 

18.52 

11.19 

0.40 

1.80 

12 
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11.19 

0.40 

1.80 

13 
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11.19 
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1.80 

14 
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0.40 

1 .80 

15 
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16 
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7.82 

0.48 

2.20 

17 

15.15 

7.82 
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15.15 

7.82 
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19! 
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23 
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35 

18.52 
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1.07 

4.01x10"® 

14,400 

0.093 

0.70 

9.34x10"^ 

42,600 

0.127 

0.43 

2.55x10"^ 

53,000 

0.102 

1.17 

1.07x10"^ 

142,000 

0.120 

0.81 

1.05x10"® 

52,00.0 

0.110 

1.27 

1.92x10"^ 

1 27 ,000 

0.170 

1.40 

. 

- 

- 

1.45 

- 

- 

0.67 

- 

- 

- 

1.64 

5.36x10"® 

10,300 

0.078 

1.11 

4.41x10"® 

17,600 

0.125 

0.41 

3.29x10"^ 

44,500 

0.068 

1.70 

3.54x10'^ 

63,700 

0.137 

0.43 

1.06x10'® 

25,000 

0.093 

2.03 

2.39x10"® 

17,600 

0.065 

0.97 

5.90x10"^ 

52 ,800 

0.059 

0.47 

2.81x10'^ 

83,200 

0.118 

0.74 

1.03x10"^ 

197,000 

0.123 

0.98 

_ 

- 

- 

1.21 

3.14x10"^ 

124,000 

0.113 

1.64 

1.27x10"^ 

260,000 

0.115 

1.09 

5.55x10"^ 

126,000 

0.137 

0.45 

_ 

- 

- 

1.56 

- 

- 

0.72 

1.05x10"^ 

593,000 

0.150 

0.48 

7.72x10"® 

160,550 

0.089 

0.52 

1.52x10"^ 

169,000 

0.100 

0.48 

Tests  25-29  are  the  same  as  Tests  11-15 
t Tests  arrested  - > 10®  cycles;  In/cycles 

NOTE:  Kp|_  = 33.33  ksl/Tn  for  all  Tests. 

2r^  “ 0.109  In  for  all  Tests. 

Overload  Preceeded  Underload  Except  in  Tests  R-12,  R-13. 


Spec. 

A-1 

A-1 

A-1 

A-2 

A-2 

A- 5 

A-3 

A- 3 

A-3 

A-10 

A-5 

A- 5 

A-6 

A-4 

A-6 

A-6 

A-10 

A-8 

A-1 3 

A-6 

A-1 3 

A- 9 

A-8 

A- 7 

A-6 

A-10 

A-8 

A-8 

A-8 

A-8 

A-13 

A-14 

A-14 

A-14 

A-1 5 

A-15 

A-1 2 

A-n 
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TABLE  II 


RESIDUAL 

STRESS 

INTENSITY 

PARAMETERS 

Test 

No. 

'^OL 

*^U0 

rMIN 

*^EFF 

.MAX 

1 

1.00 

0.11 

0.65 

+0,65 

0 

2 

1.30 

0.11 

0.42 

-0.13 

12,31 

3 

1.60 

0.11 

0.28 

-1.18 

15.54 

4 

1.90 

0.11 

0.18 

-1.40 

12,73 

5 

2.20 

0.11 

0.11 

- 

- 

6 

1.60 

-1.00 

0.22 

+0.10 

5.80 

7 

1.60 

-0,50 

0.22 

-0.26 

10.09 

8 

1.60 

0.01 

0.22 

-0.56 

12.15 

9 

1.60 

0.11 

0.22 

-0.71 

12.95 

10 

1.60 

0.22 

0.22 

-1.25 

14.84 

11 

1.80 

-1.00 

0.22 

+0.05 

6,80 

12 

1.80 

-0.50 

0.22 

+0.02 

7.16 

13 

1.80 

0,01 

0.22 

-0.78 

12.24 

14 

1.80 

0.11 

0.22 

-1.27 

13.59 

15 

1.80 

0.22 

0.22 

-1,51 

14.07 

16 

2.20 

-1,00 

0.22 

+0.01 

7,24 

17 

2.20 

-0.50 

0.22 

-0.91 

11.06 

18 

2.20 

0.01 

0.22 

- 

- 

19 

2.20 

0.11 

0.22 

- 

- 

20 

2.20 

0,22 

0.22 

- 

- 

21 

1.80 

-1,00 

0.11 

-0.42 

8.03 

22 

1.80 

-0.50 

0.11 

-0.52 

8.75 

23 

1.80 

0,01 

0.11 

-1.49 

12.54 

24 

★ 

1.80 

0.11 

0.11 

-1.47 

12.50 

30 

1.80 

-1.00 

0.30 

-0.20 

11.43 

31 

1.80 

-0.50 

0.30 

+0,21 

7.68 

32 

1.80 

0.01 

0.30 

-0.53 

12.97 

33 

1.80 

0.11 

0.30 

-0.89 

13.99 

34 

1.80 

0.22 

0.30 

-1,25 

14,73 

35 

1.80 

0.30 

0.30 

- 

- 

16-2.4 

2.40 

-1.00 

0.22 

-0,37 

9.09 

16-2.6 

2.60 

-1.00 

0.22 

-0,51 

9.19 

17-2.4 

2.40 

-0.50 

0.22 

-0,02 

7.45 

3A/D-.4a 

3.00 

-0.40 

0.10 

- 

- 

3A/D-.4b 

3.00 

-0.40 

0.10 

- 

- 

3A/D-.5 

3.00 

-0.50 

0.10 

-1.14 

7.47 

* 


Tests  25-29  are  the  same  as  Tests  11-15. 
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TABLE  III 


PREDICTED  AND  EXPERIMENTAL  RESULTS  FROM 
EXTENDED  CRACK  CLOSURE  CONCEPTS 


Test 

No. 

^UO 

EXPERIMENT 

(ksi.'^ 

PREDICTED 

(ksi/Tn) 

] 

1.00 

0.11 

0.65 

24.89 

17.76 

2 

1.30 

0.11 

0.42 

20.57 

17.10 

3 

1.60 

0.11 

0.28 

18.01 

16.43 

4 

1.90 

0.11 

0.18 

15.20 

15.76 

5i 

2.20 

o.n 

0.11 

>15.15 

15.09 

6 

1.60 

-1-00 

0.22 

14.31 

14.80 

7 

1.60 

-0.50 

0.22 

15.32 

14.60 

8 

1.60 

0.01 

0.22 

16.04 

15.79 

9 

1.60 

0.11 

0.22 

16.42 

16.43 

10 

1.60 

0.22 

0.22 

17.76 

17.36 

11 

1.80 

-1.00 

0.22 

13.22 

13.80 

12 

1.80 

-0.50 

0.22 

13.32 

13.89 

13 

1.80 

0.01 

0.22 

14.96 

15.29 

14 

1.80 

0.11 

0.22 

15.98 

15.98 

15 

1.80 

0.22 

0.22 

16.49 

16.97 

16 

2.20 

-1.00 

0.22 

11.48 

11.80 

17 

2.20 

-0.50 

0.22 

12.79 

12.39 

18t 

2.20 

0.01 

0.22 

>15.15 

14.30 

19t 

2.20 

0.11 

0.22 

H5.15 

15.09 

20t 

2.20 

0.22 

0.22 

>15.15 

16.19 

21 

1.80 

-1.00 

0.11 

12.92 

13.80 

22 

1.80 

-0.50 

0.11 

13.20 

13.89 

23 

1.80 

0.01 

0.11 

15.80 

15.29 

24 

1.80 

0.11 

0.11 

15.75 

15.98 

30 

1.80 

-1.00 

0.30 

14.84 

13.80 

31 

1.80 

-0.50 

0.30 

13.98 

13.89 

32 

1.80 

0.01 

0.30 

15.36 

15.29 

33 

1.80 

0.11 

0.30 

15.90 

15.98 

34 

1.80 

0.22 

0.30 

16.51 

16,97 

351 

1.80 

0.30 

0.30 

>18.52 

17.85 

16-2.4 

2.40 

-1.00 

0.22 

11.21 

10.80 

16-2.6 

2.60 

-1.00 

0.22 

10.69 

9.80 

17-2.4 

2.40 

-0.50 

0.22 

10.77 

11.64 

3A/D-.4at 

3.00 

-0.40 

0.10 

^11.11 

9.77 

3A/D-.4bt 

3.00 

-0.40 

0.10 

n 1 . 1 1 

9.77 

3A/D-.5 

3.00 

-0.50 

0.10 

9.09 

9.39 

* 

Tests 
+ Tests 


25-29  are  the  same  as  Tests  11-15. 
arrested  - Np  > 10^  cycles  ; 


4x10"^®  in/cycle 


dal 

da; 

^|MIN 

MIN 

’’d 

’’d 

EXPERIMENT 

PREDl 

CTED 

EXPERIMENT 

PREDICTED 

(in/cycle) 

{in/cycle) 

cycles 

cycles 

2.60x10”^ 

2.77x10"'’ 

0 

580 

3.65xl0‘^ 

2.73x10"^ 

6,800 

3,150 

3.79x10"^ 

2.11x10"® 

22,000 

20,450 

1 .86x10'^ 

6.42x10"® 

121 ,000 

262,300 

1.08x10"’® 

- 

>10® 

9.65x10'^ 

7.13x10"® 

5,850 

8,400 

5.02x10'^ 

7.86x10"® 

10.600 

7,800 

2.93x10"^ 

3.57x10"® 

11  ,400 

13,900 

2.12x10“^ 

2.11x10"® 

18,400 

20,500 

5.29x10'^ 

8.47x10"^ 

44  ,000 

34,900 

4.31x10'^ 

2.76x10"® 

1 1 ,800 

16,800 

4.01x10'^ 

2.55x10"® 

14,400 

17,800 

9. 34x10"^ 

6.36x10"^ 

42,600 

49,200 

2.55x10’^ 

2.50x10"^ 

53,000 

97,200 

1.07x10"^ 

3.75x10"® 

142,000 

388 ,800 

1 .05x10'^ 

7.37x10"^ 

52 ,000 

44,200 

1.92x10"^ 

3.46x10"^ 

127,000 

76,700 

3.67x10"^ 

- 

>10® 

-13 

1.08x10 

_ 

>10® 

. 

- 

>10® 

5.36x10"® 

2.76x10"® 

10,300 

16,800 

4.41x10"® 

2.55x10"® 

17,600 

17,800 

3.29x10"^ 

6.36x10"^ 

44,500 

49,200 

3.54x10"^ 

2.50x10"^ 

63,700 

97,200 

1.06x10'® 

2.76x10"® 

25,000 

16,800 

2.39x10'® 

2.55x10"® 

17,600 

1 7 ,800 

5.90x10'^ 

6.36x10'^ 

52,800 

49,200 

2.81x10'^ 

2.50x10"^ 

83,200 

97,200 

1.03x10'^ 

3.75x10"® 

197,000 

388,800 

_ 

1.47x10"^ 

- 

>10® 

3.14x10'^ 

5.38x10"® 

1 24 ,000 

55,600 

1 .27x10"^ 

4.93x10"® 

260,000 

59,200 

5.55x10"^ 

1.56x10"® 

126,000 

137,000 

_ 

2.15x10'® 

- 

684 ,900 

2.15x10"® 

- 

584,900 

1.05x10'^ 

5.55x10"® 

593.000 

292,000 
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time 


R 


UO 


^UL 

^OL 


*^MIN 

•^OL 


K, 


MIN 


MAX 


K 


OL 


"OL  K, 


MAX 


Figure  1.  Definition  of  variables 
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Stress  Intensity 
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Figure  2.  Q ^ vs.  IL  test  matrix 


Figure 


Stress  Intensity 


vs.  test  matrix 


Figure  8.  Typical  fracture  surfaces  for  specimens  A-5,  A-6  and  A-8. 


a (in) 


Figure  9 • Spline  function  and  integrated 

dN 

from  the  movable  strip  method  for 
a versus  N data 
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uu/u  IN  vm/oycie/ 


Figure  11.  Constant  gtress  amplitude  data  (showing  alternate 
points  only)  and  least  square  fit  line  for 
Equation  4 where  AK  = 0.5  + 0.4R 


curves  for  tests  in  the  vs.  K plane  (K^.„  = 18.52  ksi>/ta) 


w 
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cycles,  Np,  for  Q„.  vs.  R,.p,  plane 


.00  -.50  0.01  0.11  0.22  0.30 


p 
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Figure  18.  Ratio  of  measured  delay  zone  size  to  the  plane  stress 
plastic  zone  size,  Aa*/2r 


EXPERIME 


Figure  25.  Comparison  of  predic 


5 


da/dN  (in/cycle) 


